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ABSTRACT 


Y 

\  4 

The  Acoustic  Variability  Experiment  -KAVCX)  -described 
he££>  measured  acoustic  pulses  which  travel  a  wholly  refracted 
(400m)  path  between  a  source  and  receiver  in  the  upper  ocean. 
The  20  KHz  pulse  travel  times  had  rms  variations  between  1 
and  60  microseconds,  calculated  over  each  half  hour  sample 
period,  and  corresponding  pulse  amplitude  fluctuations  between 
3%  and  14%. 

1  The  temporal  and  spatial  structure  of  the  temperature 
field  were  measured  simultaneously  with  the  acoustic  trans¬ 
mission.  A  spectrum  of  a  15-day  temperature  record  showed  the 
presence  of  inertial  and  tidal  motions,  and  internal  waves 
at  higher  frequencies.  Microstructure  was  intermittent  and 
appeared  to  be  associated  with  internal  waves.  The  temperature 
integral  scales  calculated  from  the  spatial  autocorrelation 
functions  ranged  from  1.2  to  14m  (est)  in  the  vertical  and  from 
3.4m  to  greater  than  50m  (est)  in  the  horizontal.  "The  tempera¬ 
ture  structure  was  anistropic  with  the  average-hoflzontal  to 
vertical  scale  ratio  equal  to  about  fomr-rp  A  comparison  of 
the  measured  acoustic  amplitude  variance  with  temperature 
variance  and  scale  lengths  measured  at  the  receiver  station 
showed  poor  correlation  using  the  theoretical  models  of 
Chernov  and  Debye. 
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I .  INTRODUCTION 


The  Acoustic  Variability  Experiment  (AVEX)  described  here 
measured  acoustic  pulses  which  travel  a  wholly  refracted 
(400m)  path  between  source  and  receiver  in  the  upper  ocean. 

The  acoustic  pulse  was  composed  of  a  0.5ms  burst  of  20  KHz 
CW  followed  by  a  4.5ms  pseudo  random  noise  pulse  containing 
energy  in  the  frequency  range  between  5  and  25  KHz.  The 
temporal  variability  of  acoustic  amplitude  and  phase  was 
measured.  The  temporal  and  spatial  structure  of  the  tempera¬ 
ture  field  was  measured  simultaneously  with  the  acoustic 
transmission  at  the  source,  receiver  and  at  a  mid-point. 

It  is  hypothesized  that  the  acoustic  variability  for  5-25 
KHz  frequency  range  signals  in  the  upper  ocean  is  due  pri¬ 
marily  to  variations  in  the  f inestructure  of  the  index  of 
refraction  as  inferred  by  the  temperature,  and  that  the 
variations  in  the  temperature  structure  are  associated  with 
tides  and/or  internal  waves. 

There  are  a  number  of  acoustic  scattering  theories  (see 
for  example  Chernov  (1960)  and  Tatarski  (1961))  where  the 
transmission  medium  is  characterized  by  a  correlation  length 
calculated  from  the  correlation  function,  or  by  the  structure 
function,  of  the  oceanic  index  of  refraction.  The  applica¬ 
tion  of  these  theories  assumes  stationarity  and  homogeneity 
of  the  ' tatistic  describing  the  index  of  refraction.  There 
is  growing  iviuance  that  the  ocean  f inestructure  as  represented 
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by  temperature  fluctuations  is  not  homogeneous  or  stationary, 
particularly  in  the  upper  ocean  (see  for  example  Gibson,  1980) . 
It  is  expected  that  the  acoustic  variability  in  the  upper 
ocean  is  not  stationary,  but  episodic,  due  to  non-stationar- 
ity  of  the  index  of  refraction  as  measured  by  the  temperature 
structure. 

The  objectives  of  this  paper  are  to  describe  the  Acoustic 
Variability  Experiment  (AVEX) ,  and  to  perform  a  preliminary 
analysis  of  the  dependence  of  the  acoustic  phase  and  ampli¬ 
tude  fluctuations  on  the  measured  ocean  thermal  fine-  and 
microstructure  in  the  upper  ocean.  AVEX  was  conducted  in 
Carmel  Bay,  California  during  the  month  of  August  1979. 

Specific  measurements  during  AVEX  described  here  include: 

1)  Amplitude  and  phase  (travel  time)  fluctuations  of  a 
20  KHz  acoustic  pulse  transmitted  from  and  received  by 
stable,  bottom  mounted  platforms,  via  a  wholly  refracted 
path. 

2)  Temperature  f inestructure  obtained  from  horizontal 
and  vertical  thermistor  arrays  mounted  on  the  two  fixed 
bottom  platforms . 

3)  Waves  and  currents  near  the  acoustic  pulse  path,  wind 
speed  and  direction,  solar  insolation,  and  barometric  pressure 
in  the  experimental  area. 
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II.  LITERATURE  REVIEW 


Temperature  micro-  and  fine-structure  occur  in  the  ocean 
as  patches  of  thermal  inhomogeneity .  The  inhomogeneities 
are  of  various  sizes,  shapes,  and  temperatures,  which  act 
to  focus  or  defocus  the  sound  waves.  The  terms  microstruc¬ 
ture  and  finestructure  refer  to  measurement  scales  of  a  lesser 
magnitude  than  the  gross  temperature  structure.  Finestructure 
is  generally  accepted  to  reflect  vertical  structure  sizes  in 
the  order  of  meters,  while  microstructure  refers  to  fluctua¬ 
tions  with  vertical  scales  of  less  than  one  meter  (Federov, 
1978) .  It  is  necessary  to  first  separately  review  the 
microstructure  and  the  acoustics  before  the  connection 
between  the  two  can  be  investigated. 

The  stochastic  nature  of  oceanic  turbulence  must  be 
described  statistically.  Three  statistical  descriptors  to 
be  used  are:  1)  the  autocovariance  function,  2)  the  struc¬ 
ture  function,  and  3)  the  variance  spectrum.  Each  technique 
requires  a  statistically  homogeneous  medium,  and  a  certain 
degree  of  stationarity . 

1 .  Autocovariance  Function 

The  autocovariance  function,  Rx(p) ,  is  characterized 
by  the  mutual  relation  between  the  fluctuations  at  different 
times  and  positions  in  space.  It  is  expressed  by 

Rx(p)  =  E[(T(x)T(x+p) ]  (1) 
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where  E  indicates  expected  value  and  p  is  the  spatial  lag 
or  distance  between  the  two  sensors . 

2 .  Structure  Function 

The  structure  function,  Dv(p),  is  given  by 

Dx(p)  *  E[(T(x)-T(x+p))2]  (2) 

and  for  a  completely  stationary  process  reduces  to 

VP)  -  2(Rx(0)-Rx(p))  (3) 


3 .  Spectrum 

The  fourier  transform  of  the  autocovariance  function 
is  the  energy  density  spectrum  describing  the  variance  of  the 
temperature  structure  in  wave  number  or  frequency  space. 

The  one-dimensional  spectrum  in  wave  number  space  is  given 
by 


Sv(k)  *  /  Rv(p)e"27rikpdp  (4) 

X  0  * 


A.  MICROSTRUCTURE 

A  large  number  of  observations  of  thermal  structure  have 
been  made  in  recent  years  due  to  the  advent  of  improved 
instrumentation.  The  discussion  here  will  be  limited  to 
temperature  measurements  acquired  in  the  upper  ocean  per¬ 
tinent  to  AVEX.  Descriptions  of  upper  ocean  thermal  processes 
are  relatively  limited  because  the  considerable  complexity 


12 


in  this  region  of  the  ocean  does  not  lend  itself  to  a  simple 
picture . 

Observation  of  the  horizontal  and  vertical  ocean  tempera¬ 
ture  profiles  have  revealed  small  scale  variations  from  the 
gross  macroscale  changes  known  previously.  The  microscale 
inhomogeneities,  with  scales  less  than  one  meter,  exist  to 
some  extent  throughout  the  oceans,  despite  the  diffusion  and 
mixing  processes  which  tend  to  obliterate  them.  Federov 
(1978)  summarized  the  spatial  scales  of  the  f inestructure 
of  temperature  and  salinity.  He  points  out  the  strong 
anisotropy,  with  the  average  order  of  the  horizontal  to 
vertical  length  scales  on  the  order  10 . 

According  to  Munk  (1981) ,  micro  or  fines tructures  are 
due  to  internal  wave  straining  or  intrusive  processes. 
Intrusive  processes  occur  where  two  different  water  masses 
intersect  with  resulting  interleaving  as,  for  example,  occurs 
at  a  front.  The  f inestructure  is  identified  by  layering  of 
intense  temperature  inversions  which  are  generally  balanced 
by  positive  salinity  gradients,  so  that  density  increases 
stably  with  depth.  When  f inestructure  laminae  of  water  of 
differing  salinity  and  temperature  occur,  microstructure  can 
result  due  to  "fingering"  by  double  diffusion  (Turner,  1981) . 
Because  of  the  water  mass  structure  of  Carmel  Bay,  we  would 
not  normally  expect  intrusive  processes  to  be  important  in 
AVEX.  However  this  would  not  be  true  in  winter  when  the 
Carmel  River  discharges  fresh  water  into  the  Bay. 
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The  more  commonly  observed  (Woods,  1968)  or  inferred 
mechanism  for  microstructure  generation  is  due  to  internal 
waves.  The  occurrence  of  internal  wave  breaking  has  been 
characterized  by  Munk  (1981)  as  analogous  to  the  intermittent 
white  capping  of  surface  waves  in  a  light  wind.  The  inter¬ 
mittent,  or  episodic  nature,  of  breaking  internal  waves  re¬ 
sults  in  the  observed  patchiness  of  microstructure. 

Because  internal  wave  breaking  is  episodic,  the  rela¬ 
tively  high  frequency  induced  temperature  fluctuations  are 
nonstationary  over  time  (length)  scales  of  the  order  of 
several  internal  wave  periods  (lengths) . 

Munk  (1981)  gives  an  interesting  summary  of  the  present 
knowledge  of  the  little  understood  subject  of  breaking  internal 
waves.  Breaking  can  apparently  be  due  to  either  advective 
instability  or  shear  instability.  Advective  instability 
occurs  when  the  water  particle  speed  at  the  crest  exceeds 
the  wave  speed,  similar  to  wave  breaking  on  a  beach.  The 
breaking  condition  occurs  for  large  amplitude,  or  large 
slope,  internal  waves,  and  can  occur  in  the  absence  of  shear, 
although  it  is  augmented  by  shear. 

The  limiting  case  of  shear  instability  is  Kelvin-Helmholtz 
instability  which  can  occur  at  an  abrupt  density  transition 
even  in  the  absence  of  finite  amplitude  waves.  Munk  (1981) 
provides  an  argument  that,  for  linear  internal  waves  only 
in  the  presence  of  self-shear,  conditions  favor  advective 
instability  over  shear  instability,  and  that  the  wave  field 
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is  within  a  small  numerical  factor  of  advective  insta¬ 
bility.  The  argument  for  shear  instability  can  be  increased 
by  considering  a  spectrum  of  internal  waves  where  the  lower 
frequency  modes  impose  an  ambient  shear  environment  for  the 
higher  modes  that  enhances  breaking  conditions.  Taking  into 
account  a  spectrum  of  internal  waves,  it  is  suggested  that 
instabilities  can  occur  due  either  to  advection  or  shear, 
depending  on  wave  slope  and  the  ambient  shear. 

The  wave  number  spectrum  of  horizontal  temperature  varia¬ 
tions  can  be  divided  into  several  ranges  characterized  by 
spectral  slope  power  laws.  Garrett  and  Munk  (1972,  1975, 

1979)  proposed  evolving  semi-empirical  spectral  models  that 
have  provided  impetus  for  other  investigators  to  acquire 
data  to  compare  with  the  models.  The  Garrett  and  Munk  theory 
predicts  a  -2  power  law  decay  at  lower  frequencies  and  a 
-2.5  law  at  higher  frequencies.  The  internal  wave  contribution 
cuts  off  at  the  local  Vaisala  frequency,  N.  Above  N,  Levine 
and  Irish  (1981)  have  shown  a  continuation  of  the  -2.5  slope 
attributed  to  f inestructure  alone.  At  larger  wave  numbers, 
buoyancy  and  convective  ranges  are  found  with  -3  and  -5/3 
power  laws  respectively. 

Kolmogorov  (1941)  theorized  that  energy  was  extracted 
from  an  unstable  mean  flow  by  the  largest  eddies  in  a  turbu¬ 
lent  flow.  The  turbulent  flow,  consisting  of  eddies  of 
various  sizes,  progressively  transferred  the  energy  to  the 
smaller  eddies  where  heat  was  ultimately  dissipated  by 
action  of  viscosity  in  the  smallest  eddies.  From  this  theory. 
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Kolmogorov  predicted  a  2/3  power  law  dependence  for  the 
structure  function  for  turbulent  flow.  The  transformation 
of  this  result  to  wave  number  space  energy  spectra  yields 
a  wave  number  to  the  -5/3  power  law  dependence.  This  was 
first  demonstrated  by  Batchelor  (1953)  and  later  verified 
by  Grant  et  al.  (1962)  in  the  ocean  at  high  wave  numbers. 

In  recent  years,  there  has  been  a  veritable  explosion 
in  the  number  of  fine-scale  temperature  measurements  due 
to  the  advent  of  improved  instrumentation.  Most  of  the 
measurements  use  profiling  devices  and  are  concerned  with 
vertical  scales.  Federov  (1973)  gives  a  good  summary  up 
to  1974,  including  much  of  the  Russian  literature.  More 
recent  reviews  and  references  include  Phillips  (1977) , 

Turner  (1981) ,  Munk  (1981) ,  Gargett  and  Osborn  (1981) . 
Greatest  attention  has  been  given  to  comparing  the  tempera¬ 
ture  variability  about  the  mean  profile  and  comparing  the 
results  with  the  Garrett-Munk  spectra.  Due  to  stationarity 
constraint,  the  upper  ocean  is  usually  not  included  in 
the  spectral  calculations. 

Horizontal  microstructure  scales  in  the  upper  ocean  have 
been  measured  by  thermistors  mounted  on  submarines  and  towed 
fish.  As  early  as  1948,  Urich  and  Searfoss  utilized  a  sub¬ 
marine  to  measure  microstructure  off  the  Florida  coast  at 
depths  up  to  50m,  and  they  measured  length  scales  of  20  to 
100cm  and  typical  temperature  variations  of  0.02  to  0.1° C. 
Leiberman  (1951)  made  similar  measurements  off  Southern 
California  in  the  depth  range  of  30  to  60m  and  calculated 
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mean  patch  size  on  the  order  of  60cm  and  an  average  tempera¬ 
ture  variation  of  0.05°C.  Zenk  and  Katz  (1975)  used  a  towed 
thermistor  c.h  26m  depth  in  the  Sargosso  Sea.  They  point  out 
that  the  use  of  spectral  analysis  in  a  high-frequency  range 
may  be  seriously  affected  by  the  nonstationarity  of  the 
horizontal  temperature  data  series  due  to  the  patchiness  of 
the  temperature  variations  encountered.  Zenk  and  Katz 
(1975)  calculated  37  sequential  spectra,  each  spectrum  repre¬ 
senting  approximately  lKm.  They  found  the  spectral  density 
varied  by  greater  than  30  fold,  which  is  consistent  with 
earlier  Russian  work  reported  in  Zenk  and  Katz  (1975)  .  More 
recently,  Gibson  (1980)  measured  both  temperature  and  velocity 
variations  with  a  towed  fish.  He,  and  also,  Crawford  and 
Osborn  (1980)  using  a  free-fall  profiler  measuring  both 
temperature  and  shear  velocity,  found  most  often  that  small 
scale  temperature  inversions  and  velocity  microstructure 
occur  together.  Less  frequently,  Gibson  (1980)  found  patches 
of  temperature  microstructure  without  velocity  microstructure. 
These  patches  appear  to  be  the  remains  of  internal  waves  that 
had  previously  broken  and  are  referred  to  as  "fossil  turbulence". 

Munk  (1981)  writes,  "the  picture  that  emerges  is  one  of 
a  fines true ture  that  is  dominated  by  internal  wave  straining 
and  is  fairly  uniform,  in  contrast  to  microstructure  that 
is  extremely  patchy  and  variable  even  in  the  mean." 

B.  ACOUSTIC  PROPAGATION 

To  describe  wave  propagation  through  the  ocean,  the  basic 
Helmholtz  wave  equation  for  homogeneous  media  must  be  modified 
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to  account  for  changes  of  the  refractive  index  with  space  and 
time.  Historically,  much  of  the  research  into  wave  propaga¬ 
tion  through  fluctuating  media  has  been  derived  from  theoreti¬ 
cal  studies  of  electromagnetic  wave  propagation  in  the 
atmosphere.  Treatments  based  on  work  by  Rytov  (1937) , 

Chernov  (1960)  and  Mintzer  (1953)  have  included  refractive 
index  fluctuations  into  the  Helmholtz  wave  equation  under 
different  assumptions  of  isotropy,  scale-size,  and  the  strength 
of  inhomogeneities  within  the  media.  The  methods  and  predic¬ 
tions  of  acoustic  amplitude  variations  from  this  early  work 
are  described  below,  together  with  a  summary  of  recent  exten¬ 
sions  of  propagation  theory  including  de  Wolf  (1975)  and 
Flatte  (1978). 

Mintzer  (1953)  applied  the  Born  approximation  to  the 
basic  wave  equation.  The  Born  approximation  is  a  perturba¬ 
tion  technique  which  is  valid  only  if  the  acoustic  amplitude 
fluctuations  remain  small,  and  accounts  for  scattering  only 
from  a  single  turbulent  blob  within  the  acoustic  path  length. 
The  solution  yields  an  expression  relating  the  form  of  the 
structure  function  of  the  scattering  medium  to  the  coefficient 
of  variation  of  acoustic  pulse  amplitude,  V.  An  assumption 
used  in  the  following  three  treatments  is  that  during  the 
travel  time  of  an  acoustic  pulse,  the  medium  remains  constant 
(i.e.,  the  microstructure  mozaic  remains  fixed)  and  changes 
in  the  microstructure  occur  only  between  successive  pulses. 

The  following  equations  and  conditions  were  presented  by 
Stone  and  Mintzer  (1962,  1965)  after  they  conducted  a  series 
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of  laboratory  experiments  to  check  predictions  of  the  vari¬ 
ance  of  normalized  sound  pressure  level  in  a  low  frequency 
(wave  equation  vs.  ray  theory)  region: 


V2  ■  j  /it  Kq  a2  a  L  (Gaussian  rx(p))  (5) 

2  2  2 

V  =  oi  Kq  a  L  (Exponential  R  (p))  (6) 

where 

V2  «  E[(^)2]  (7) 

p0 

p  =  Acoustic  pressure 

Pq  =  mean  acoustic  pressure 

a  =  Taylor  structure  length  scale 

=  1.14  x  Integral  length  scale  (for 

Gaussian  correlation  function) 

L  -  Acoustic  path  length 

2 

a  =  Variance  of  index  of  refraction 

A  =  Acoustic  wavelength 

Kq  =  2tt/\  (8) 

and 

A  <<  27ra  (9) 


A  <<'  2ttL 


(10) 
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The  wave  region  applies  when 


a  >>  2a2/L  (11) 

Mintzer  applied  Liebermann' s  (1951)  experimental  values  to 
his  results  and  also  obtained  good  agreement  between  the 
theoretical  results  and  ocean  measurements  obtained  by 
Sheehy  (1950) .  The  severe  wavelength  restrictions  required 
by  this  theory  limited  its  general  use. 

Chernov  (1960)  developed  formulae  relating  acoustic 
amplitude  fluctuations  to  path  length  and  microstructure 
scale  lengths  for  both  ray  theory  and  wave  theory  regimes. 

The  ray  theory,  or  geometrical  optics  approach,  is  a  lineari¬ 
zation  of  the  wave  equation,  and  employs  the  assumptions  that 

the  scale  of  the  inhomogeneities  is  large  compared  to  the 

2 

acoustic  wavelength,  \  <<  £q/L,  where  £q  is  the  inner  scale 
of  the  turbulence  and  transit  time  is  small  compared  to  the 
rate  of  change  in  the  refractive  index  field.  This  approach 
ignores  diffractive  effects  and  only  considers  Rayleigh 
scattering  and  refractive  effects.  Liebermann  (1951) ,  apply¬ 
ing  this  theory,  predicted  variations  in  sound  amplitude 
with  the  following  equations: 

vi  =  A  *2  (L/a)3  (12) 


where 
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(13) 


V, 


E“f» 


and 


E[(^)2) 


(14) 


Amplitude  variations  measured  by  Liebermann  were  larger  than 
those  predicted  due  to  temperature  inhomogeneities.  Lieber¬ 
mann  suggested  that  additional  contributions  to  the  fluctua¬ 
tions  could  be  due  to  biological  scattering,  chemical 
relaxation  and  bubbles  in  the  near-surface  region.  Upper 
ocean  acoustic  measurements  made  by  Sagar  (1955,  1957,  1960) 
showed  that  sea  state  was  also  a  contributing  element  to 
amplitude  fluctuations. 

In  the  wave  theory  regime,  the  method  of  small  perturba¬ 
tions  developed  by  Chernov  requires  that  the  amplitude  of 
the  scattered  waves  is  "Small  compared  to  the  amplitude  of 
the  initial  wave,  thereby  limiting  the  range  in  which  this 
method  can  be  applied.  Since  large  amplitude  and  phase  fluc¬ 
tuations  are  often  observed  in  situations  where  wave  theory 
is  still  applicable,  an  alternate  method  (Rytov)  was  employed 
by  Chernov  to  remove  the  range  restriction. 

Rytov' s  method  of  smooth  perturbations  was  used  by 
Tatarski  (1961)  to  investigate  fluctuations  of  electromagnetic 
and  acoustic  waves  by  developing  a  stochastic  wave  model  for 
the  atmosphere.  It  includes  the  geometrical  optics  and  Born 
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approximation  method  as  special  cases.  The  restrictions  of 
the  model  are  that  the  scale  of  turbulence,  % ,  is  much  larger 
than  the  acoustic  wavelength  and  that  the  medium  is  weakly 
inhomogeneous  and  quasi-stationary . 

The  Debye  approximation  to  the  wave  equation  is  a  general 
model  which  describes  the  turbulent  ocean  effects  on  sound 
transmission  without  special  long  range  or  high-frequency 
restrictions  (Neubert  1970).  It  contains  Rytov's  method  as 
part  of  its  solution  and  reduces  to  the  Born  (single  scatter¬ 
ing)  model  for 


A  >>  2  it  ot/aL  (15) 

It  predicts  a  variance 

V2  =  2  a2  K2  l  L  (16) 

where  l  is  the  turbulent  integral  scale. 

Much  of  the  recent  work  in  acoustic  propagation  in  the 
ocean  has  concentrated  on  long  range  propagation  experiments , 
in  which  the  use  of  acoustics  to  remotely  probe  the  ocean  has 
been  investigated.  This  research  prompted  the  development 
of  a  more  unified  treatment  of  both  weak  and  strong  scatter¬ 
ing  regimes.  The  different  scattering  regimes  for  acoustic 
transmission  in  the  ocean  have  been  mapped  into  a  strength/ 
diffraction  parameter  space  by  de  Wolf  (1975)  and  Flatte  (1979) . 
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The  boundaries  within  which  the  different  propagation  theories 
apply  may  be  defined  in  terms  of  strength  and  size  parameters 
which  characterize  the  inhomogeneities  in  the  ocean.  For 
the  isotropic  case,  Flatte  (1979,  p.  91)  has  defined  the 
dimensionless  size  parameter  A  as 


A  =  L_1  /  (2tt)  _1  [R?  (x)  /L2]  dx 


Evaluating  the  integral 


A  =  L/(6  £  Kq) 


where  R^(x)  is  the  radius  of  the  first  Fresnel  zone  for  a 
source  and  receiver  separated  by  a  range  L,  and  l  is  the 
correlation  length  of  the  temperature  fluctuations  and  x  is 
the  range  from  the  source. 

The  ensemble  average  integral  of  the  sound  speed  along  a 
straight  line  between  the  source  and  receiver  is  used  to 
characterize  the  strength: 


$2  -  <(Kn  /  u  dx) 2> 


where 


=  the  strength  parameter  (dimensionless) 

=  the  acoustic  wavelength 

*  the  fractional  sound-speed  deviation 
from  the  mean  sound  channel  value 
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For  the  case  where  >>  l,  $  can  be  approximately  evaluated  as 

2  2  2 

$  =  qQ  <y  >  L  (20) 

where  z  0.4  is  the  integral  scale  length.  Under  the 

2 

geometrical  optics  approximation  the  quantity  $  is  the  rms 
variation  in  phase  of  a  pulse  at  the  receiver. 

Several  important  regime  boundaries  are  shown  in  the 
representation  of  A-$  space  shown  in  Fig.  1  (after  Flatte 
1979,  where  a  full  discussion  may  be  found).  The  geometri¬ 
cal  optics  and  Rytov  extension  regions  have  been  discussed 
earlier.  The  most  important  boundary  in  the  diffraction 
region  (A  >  1)  is  the  line  4>  »  1,  where  the  rms  phase 
fluctuations  and  the  intensity  fluctuations  both  approach 
unity.  Here  phase  fluctuations  are  normalized  by  the 
acoustic  wave  period  (50us  for  20KHz) ,  and  acoustic  inten¬ 
sity  is  normalized  by  its  mean  value. 

The  measurements  made  in  the  AVEX  experiment  will  be 
discussed  in  terms  of  their  position  in  this  parameter  space. 

C.  ACOUSTIC  PROPAGATION  EXPERIMENTS 

During  the  1950's  several  ocean  measurements  of  acoustic 
amplitude  were  made  from  drifting  ships;  see  for  example 
Sheehy  (1950),  Liebermann  (1951)  and  Sagar  (1955  and  1957). 
These  authors  found  amplitude  fluctuations  significantly 
larger  than  their  theoretical  predictions  and  suggested  that 
this  was  due  to  environmental  effects  other  than  temperature 
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Figure  1.  Regimes  of  a-*  space 
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inhomogenf;ities ,  for  example  bioscatterers  and  air  bubbles, 
and  measurement  problems  such  as  rotation  of  non-omnidirectional 
transducers  and  interference  from  surface  and  bottom  reflected 
pulse  arrivals.  The  use  of  drifting  platforms  did,  however, 
allow  estimates  of  range  dependence  to  be  made  relatively 
easily. 

Examples  of  recent  long  range  acoustic  experiments  include 
the  MIMI  project  (see  for  example  Steinberg,  et  al.  (1972)  and 
Clark  and  Weinberg  (1973))  during  which  406  Hz  pulses  were 
transmitted  over  a  1250Km  wholly  refracted  path.  This  long 
range  placed  the  experiment  within  the  saturated  area  of 
Flatte's  parameter  space;  i.e.,  there  were  several  rays  or 
micro-multipaths  from  the  source  to  the  receiver. 

Porter  and  Spindel  (1977)  report  on  a  propagation  experi¬ 
ment  using  a  ship  suspended  source  and  free  drifting  receiver 
buoys.  An  acoustic  doppler  navigation  system  was  used  to 
measure  buoy  motion.  They  calculated  frequency  spectra  of 
acoustic  phase  for  hydrophones  at  depths  between  300  and  1500m 
and  found  that  the  shallowest  receivers  exhibited  an  order  of 
magnitude  greater  phase  fluctuations  than  predicted  by  a 
random  multipath  model  based  on  sound  interaction  with  internal 
waves. 

Kennedy  (1969)  measured  arrival  times  and  amplitudes  of 
10ms,  800Hz  pulses  transmitted  over  a  25-mile,  wholly 
refracted  path.  Fixed  sources  were  utilized  with  an  acrustiv. 
array  of  11  vertically  spaced  hydrophones  with  a  maximum 
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separation  of  300  feet.  Wherever  possible,  the  measured  re¬ 
sults  were  compared  with  the  theoretical  values  of  Chernov 
(1960) .  Results  calculated  for  phase  and  amplitude  fluctua¬ 
tions  compared  favorably.  Kennedy  found  evidence  that  some 
of  the  amplitude  fluctuations  were  caused  by  the  interference 
of  secondary  scattered  waves  with  the  primary  wave.  Calcu¬ 
lations  were  made  using  the  theory  of  an  inhomogeneous  medium 
with  multiple  patch  sizes  as  described  by  Kolmogorov  (1941)  . 

Ewart  (1976)  measured  amplitude  and  phase  fluctuations 
in  pulses  sent  over  a  wholly  refracted  path  over  a  range 
17.2Km  and  depth  of  1000  meters.  Eight-cycle  pulses  at 
4, 166Hz  and  16-cycle  pulses  at  8,333Hz  were  sent  alternately 
every  15.7  seconds  and  received  by  three  horizontally  spaced 
hydrophones.  Power  spectra  of  amplitude  and  phase  fluctua¬ 
tions  at  single  receivers  and  phase  differences  for  the 
series  of  three  receiver  hydrophones  were  calculated  and 
compared  with  theoretical  calculations .  The  amplitude 
fluctuation  spectra  displayed  strong  diurnal  and  semidiurnal 
peaks,  while  the  phase  fluctuation  spectra  displayed  a  third 
peak  at  the  quarter  diurnal  frequency.  In  the  spectral  band 
between  the  inertial  frequency  for  the  local  area  and  the 
Vaisala  frequency  at  1000  meters,  Ewart  found  power  law 
dependencies  with  frequency,  ut,  of  oo  \  and  for 

phase,  amplitude  and  phase  difference  spectra,  respectively. 

In  a  subsequent  analysis  of  Ewarts  data,  Oesaubies  (1976) 
attributed  the  source  of  phase  fluctuations  with  frequencies 
between  the  inertial  and  buoyancy  frequencies  to  be  caused 
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by  internal  wave  activity.  He  was  able  to  use  the  geometric 
approximation  in  this  model,  and  emphasized  that  the  model 
was  very  sensitive  to  oceanic  measurements. 

The  contribution  of  fine  to  micro-scale  temperature 
activity  to  acoustic  amplitude  and  phase  fluctuations  was 
not  well  defined  by  the  experiments  described  above.  Desaubies 
(1976)  did  point  out  that  phase  fluctuations  are  dominated 
by  larger  scale  sound-velocity  variations,  whereas  the 
amplitude  fluctuations  are  related  to  smaller  scales.  He 
postulated  that  the  amplitude  variations  are  not  caused 
directly  by  the  internal  wave  field,  but  by  small  scale  multi¬ 
paths  caused  by  the  temperature  f inestructure  being  advected 
by  the  internal-wave  field. 
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III.  ACOUSTIC  VARIABILITY  EXPERIMENT  (AVEX) 


A.  EXPERIMENT  DESCRIPTION 

Acoustic  and  oceanic  data  were  acquired  from  August  2 
to  August  26,  1979  in  the  Carmel  Canyon  of  Carmel  Bay.  Figure 

2  is  a  schematic  cross-section  of  the  experiment,  and  Figure 

3  shows  a  plan  view.  The  bathymetry  along  the  acoustic  path 
had  relatively  level  terrain  at  depths  between  30  and  35m  each 
side  of  the  Canyon  and  steep  slopes  to  the  bottom  depth  of 
180m.  A  fathometer  trace  of  the  bottom  profile  from  the 
transmitter  station  to  hydrophone  station  400m  away  is  shown 
in  Figure  4.  The  location  was  chosen  because  the  bathymetry 
offered  the  opportunity  to  have  a  bottom-mounted,  transmitter- 
transducer  and  receiver  hydrophones  located  on  each  side  of 
the  Canyon,  with  no  surface  or  bottom  reflected  multipath 
arrivals  for  the  duration  of  the  signal  pulse.  The  use  of 
the  rigid,  bottom-mounted  transducer  and  hydrophones  precluded 
any  biasing  of  the  amplitude  and  phase  fluctuations  due  to 
motion  of  the  transducer  or  hydrophones .  The  experiment  can 
be  segmented  into  five  sites:  1)  Station  T  (Transmitter); 

2)  Station  R  (Receiver) ;  3)  Station  P  (Profiler) ;  4)  Current 
Meter  Array;  and  5)  Shore  Station.  The  five  stations  are 
described  separately. 

1.  Station  T 

The  transmitter  tower  is  shown  schematically  in  Figure 
5,  and  Figure  6  is  a  photograph  of  the  deployment  of  the 
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Figure  2.  Schematic  cross-section  of  experiment  location. 
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Figure  4.  Fathometer  trace  across  Carmel  Canyon. 
Depth  in  fathoms. 
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Figure  6.  Photograph  of  transmitter  tower. 


tower.  A  spherical,  omnidirectional  transmitter-transducer 
was  mounted  at  the  top  of  the  rigid  tripod,  7m  from  the  base. 
An  omnidirectional  transducer  was  used  to  ensure  the  trans¬ 
mission  of  broad-band  signals  without  the  formation  of  narrow, 
directional  acoustic  beams,  which  would  otherwise  make  the 
received  signal  amplitude  very  susceptible  to  very  small 
transducer  rotations.  The  tripod  was  located  in  32m  of  water 
on  relatively  level  terrain.  A  thermistor  array  was  attached 
to  the  legs  of  the  tripod  4.2m  above  the  bottom.  The  array 
consisted  of  eight  horizontal  and  two  vertical  thermistors 
with  one  thermistor  in  common.  The  thermistors  were  unevenly 
spaced  along  the  arrays  to  maximize  the  number  of  spatial  lags 
for  computing  correlation  and  structure  functions .  Table  I 
lists  the  spacings  for  the  three  thermistor  arrays  used  in 
the  experiment.  Two  Bendix  current  meters  were  positioned 
on  one  cross-member  of  the  tripod  to  measure  orthogonal  com¬ 
ponents  of  the  current.  A  100  psig  pressure  transducer  was 
also  mounted  on  the  tripod  to  measure  tide  and  long  period 
wave  formation.  The  oceanic  signals  measured  at  Station  T 
were  PCM  (Pulse  Code  Modulation)  encoded  and  cabled  back 
to  shore  for  recording,  using  half  of  the  4-conductor  cable 
which  also  carried  the  lKw  signal  pulses  to  the  acoustic 
transducer. 

2 .  Station  R 

Figure  7  is  a  schematic  view  of  the  Receiver  Station 
tripod,  including  the  thermistor  array.  Figure  8  is  a  photo¬ 
graph  taken  during  installation.  The  receiver  tripod  was 
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Figure  7.  Receiver  tower  schematic. 
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Figure  8.  Photograph  of  receiver  tower. 


Table  I 


Thermistor  Array  Spacings  (cm) 
Thermistor  Spacings  (cm) 


Element 

1  2 

3 

4 

5 

6 

7 

8 

9 

10 

11  12 

Station  T 

48.5 

125 

40 

65 

5.5 

10.5 

20 

84.5 

Station  R 

125.5 

5 

10 

61.5 

85 

19.5 

10.5 

5.5 

65 

40  125.5 

Station  P 

5 

60 

50 

30 

70 

10 

5 

20 

35 

10  55 

located  in  35m  of  water  on  a  slope  of  approximately  10  degrees, 
which  required  minor  adjustments  by  divers  to  level  the  horizon¬ 
tal  thermistor  array.  Eight  thermistors  were  used  in  the 
horizontal  array,  and  five  in  the  vertical  array,  with  one 
common  to  both.  The  two-hydrophone  array  mounted  on  top  of 
the  tripod  was  7m  above  the  bottom.  The  hydrophones  were 
spaced  4.15m  apart  in  the  horizontal  and  oriented  approxi¬ 
mately  perpendicular  to  the  transmission  path.  Acoustic 
signals  received  at  the  station  were  VHF ,  FM  telemetered  to 
shore.  The  measured  oceanic  data  (wave  height  and  temperature) 
were  PCM  encoded  and  also  telemetered  to  the  shore . 

A  36m  vertical  spar  buoy  located  15m  south  of  the 
tripod  supported  the  PCM  telemetry,  four-broadband  hydrophone 
telemetry  transmitters,  two  VHF  antennas,  a  two-hydrophone 
array  and  a  wave  gage.  The  spar  buoy  was  attached  to  a  bot¬ 
tom  anchor  by  a  truck  universal  joint  which  allowed  the  spar 
to  free-swivel  (but  not  rotate) .  The  mid-depth,  two- 
hydrophone,  horizontal  array  was  located  on  the  spar  buoy, 

13.7m  below  the  mean  water  surface  with  a  spacing  of  3.86m. 
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A  pressure  transducer  used  to  measure  waves  was  positioned 
just  below  the  lower-low  water  mark  on  the  spar  buoy.  The 
station  was  powered  by  a  220  AH  18  VDC,  lead-acid  battery 
pack  mounted  on  the  tripod. 

3 .  Station  P 

A  yo-yo  profiler  was  located  near  the  middle  of  the 
Canyon  approximately  100m  east  of  the  direct  acoustic  path, 
as  shown  schematically  in  Figure  2.  The  profiler  was  de¬ 
signed  to  make  repeated,  high-resolution  density  profiles 
using  a  horizontal  array  of  12  thermistors,  a  pressure  sensor 
and  conductivity  cell,  by  moving  up  and  down  a  taut-wire 
mooring  cable  between  the  depths  of  10  and  100m  at  approxi¬ 
mately  20  cm/sec,  once  every  two  hours.  Profiling  was 
accomplished  by  varying  the  buoyancy  by  -2  percent  of  neutral, 
utilizing  a  timed  variable  ballast  device.  The  oceanic  sig¬ 
nals  were  PCM  encoded,  then  inductively  coupled  to  the  mooring 
cable.  A  surface  PCM  telemetry  buoy  was  tethered  from  the 
subsurface  profiler  buoy  and  used  to  VHF  telemeter  the  PCM 
signal  from  the  mooring  cable  to  the  shore  station.  The 
yo-yo  profiler  was  damaged  during  installation  and  did  not 
function  during  the  experiment.  The  lack  of  thermal  micro¬ 
structure  at  midrange  consequently  limited  the  results. 

4 .  Current  Meter  Array 

A  two-Aanderra  current  meter  array  was  located  approxi¬ 
mately  30m  south  of  the  profiler,  anchored  at  a  depth  of 
100m  (see  Figures  1  &  2) .  The  current  meters  were  15  and  46m 
below  the  mean  surface,  bracketing  the  sound  path.  Current 
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speed  and  direction  data  were  recorded  internally  in  digital 
format  at  5-minute  intervals  for  the  26-day  period.  The 
array  was  retrieved  by  the  use  of  timed  releases. 

5 .  Shore  Station 

A  mobile  bus,  serving  as  the  shore  station,  housed 
all  the  receiving,  pre-processing  and  recording  equipment. 
Figure  9  is  a  photograph  of  the  10m  weather  tower  attached 
to  the  bus.  The  environmental  instruments  measured  insola¬ 
tion,  back-scattered  radiation,  wind  speed  and  direction, 
and  barometric  pressure.  Signals  were  received  at  the  bus 
via  PCM  or  FM  telemetry  from  Stations  P  and  R,  and  via  the 
armored  cable  from  Station  T. 

B.  INSTRUMENTATION 

Table  II  is  a  summary  of  the  measurement  devices  utilized 
in  the  experiment  and  pertinent  information  on  each  instru¬ 
ment. 

Sound  pulses  were  emitted  by  an  omnidirectional  trans¬ 
ducer  type  C7JP3.  Omnidirectional  hydrophones  from  AN/SSQ57A 
sonotuoys  were  used  as  receivers . 

Temperature  was  measured  using  Victory  Engineering  Corpora¬ 
tion  .020  inch  bead-m-glass  probe  thermistors.  The  thermis¬ 
tors  had  a  still  water  time  constant  of  55ms  and  a  dissipation 
constant  of  0.75  MW/8C.  A  Hewlett-Packard  HP2801A  quartz 
thermometer  with  an  accuracy  of  t  .005°C  was  used  to  calibrate 
the  thermistors  through  their  associated  preconditioning 
electronics  and  PCM  encoder.  Calibration  was  accomplished 
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Figure  9.  Photograph  of  shore  station. 


Table  II.  Instrument  Summary 
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(shore  station)  search,  Inc. 

Pressure  Transducer  Transducer  Inc.  5AP-69F  0-50  psia  Nonlinearity  S 

(spar  buoy)  hysterysis 


by  cooling  a  3  liter  volume  of  water  down  to  3°C.  Thermistors 
from  each  array  together  with  the  reference  quartz  thermometer 
were  placed  in  the  stirred  volume  of  water  and  the  system  was 
allowed  to  slowly  warm  to  20  °C  (ambient  temperature) .  The 
thermistors  were  connected  to  an  HP-9845  computer  via  '  neir 
preconditioning  electronics,  PCM  encoder,  and  PCM  decoder. 

The  computer  recorded  and  stored  measurements  every  0.5°C. 

A  least-squares  linear  fit  over  the  temperature  range  of  8-16 
°C  was  then  calculated  for  each  thermistor. 

Currents  were  measured  by  both  ducted  impeller  and 
Savonious  current  meters.  The  ducted  impeller  current  meters, 
mounted  on  the  transmitter  tower,  were  Bendix  Corporation, 
model  B-10 .  The  impeller  assembly  had  a  starting  threshold 
of  2  cm/s,  a  threshold  of  linearity  of  3  cm/s,  and  a  response 
time  less  than  0.1  second.  The  impeller  rotations  were 
sensed  through  a  magnetic  read  switch.  Flat  impeller  blades 
were  used  for  equal  responses  for  both  forward  and  backward 
flow.  Aanderaa  meters,  model  RCM4 ,  were  used  to  measure  cur¬ 
rents  at  the  mid-Canyon  current  meter  array.  These  recording 
meters  use  a  Savonious  rotor  current  velocity  sensor  and  a 
magnetic  compass  for  direction  determination. 

Sea  surface  elevation  was  measured  using  a  Transducer, 
Inc.,  pressure  transducer,  model  5AP-69F,  mounted  on  the  re¬ 
ceiver  station  spar  buoy. 

Insolation  and  back-radiation  were  measured  using 
two,  model  PSP,  Eppley  precision  pyranometers .  Wind  speed 
and  direction  were  measured  using  a  wind  speed  sensor,  model 
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1022S,  and  a  wind  direction  sensor,  model  1022D,  manufactured 
by  Meteorology  Research,  Inc. 

C.  PROCEDURE 

Acoustic  and  oceanic  data  were  collected  for  28  min  every 
two  hours  starting  7  PM  August  2  to  August  26,  1979.  The 
acoustic  pulses  were  transmitted  from  Station  T  transducer, 
received  at  the  Station  R  hydrophones,  and  telemetered  to 
shore  where  they  were  processed  and  recorded.  The  measure¬ 
ment  cycle  consisted  of  the  transmission  and  reception  of 
ten  6  KHz  sinusoidal  pulses  immediately  before  and  after  a 
series  of  300  composite  pulses,  each  at  a  repetition  rate  of 
5  seconds.  Complex  spectra  of  three,  one-minute  average, 

4-25  KHz  ambient  noise  spectra  received  by  hydrophones  HI  and 
H2  were  recorded  at  the  end  of  the  pulse  measurements. 

The  6  KHz  acoustic  pulses  were  used  to  assist  in  the  de¬ 
tection  of  multipathing.  The  300  composite  pulses  were 
analyzed  to  obtain  phase  and  amplitude  fluctuations  of  the 
20  KHz  sine  wave  and  broadband  signal  contained  in  each  pulse. 
The  5-second  pulse  repetition  rate  was  determined  principally 
by  the  data  transfer  and  processing  times,  and  was  long  enough 
for  the  reflected  and  reverberated  signals  of  one  pulse  to 
decay  before  the  transmission  of  the  next  pulse. 

An  example  of  the  composite  pulse  is  shown  in  Figure  10. 
The  first  0.6  ms  of  the  pulse  consisted  of  12  cycles  of  a  20 
KHz,  zero  phase  start  and  zero  phase  finish,  sinusoidal  wave, 
followed  by  4.5  ms  of  pseudo-random  noise.  The  pseudo-random 
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Figure  10.  Current  wave  form  of  the  composite  signal 
pulse  (upper  panel),  and  voltage  wave  form  of  the  com 
oosite  signal  pulse  (lower  panel). 


signal  had  a  flat  frequency  spectrum  from  5  to  25  KHz.  The 
signal  was  repetitious  and  identical  for  each  composite 
pulse.  Each  time  the  HP-9845  computer  activated  the  pulse 
generator,  the  digitizing  delay  generator  and  timer  units 
were  initiated.  The  signal  was  fed  from  the  pulse  generator 
to  a  1  Kw,  300  VRMS  power  amplifier.  The  amplifier  output 
provided  a  5  to  25  KHz  bandwidth  constant  current  source  to 
the  transducer  via  the  shore  cable.  The  high  powered  ampli¬ 
fier  was  required  because  of  the  need  to  use  an  omnidirec¬ 
tional  transducer  while  maintaining  an  acceptable  signal  to 
noise  ratio  at  the  receiver  hydrophones.  Frequencies  below 
the  5-25  KHz  pseudo  random  noise  band  were  filtered  out  to 
prevent  clipping  in  the  constant-current  amplifier.  The  1  KW 
output  amplifier  was  used  in  a  constant  current  mode  to  enable 
the  transmitted  spectrum  to  equal  the  input  spectrum,  despite  the 

complex  impedance  characteristic  of  the  cable  and  transducer. 

The  "acoustic  pulses  were  received  by  the  four  hydrophones 
at  Station  R.  The  receiver  station  was  synchronized  to  the 
transmitter  so  that  only  the  5  ms  sound  pulse  would  be  re¬ 
corded.  Synchronization  was  accomplished  by  establishing  the 
pulse  generation  times  and  delaying  the  recording  by  the  256 
ms  required  for  the  pulse  to  travel  the  wholly  refracted  path. 

The  received  signals  were  first  fed  into  underwater  preamp 
units  and  then  telemetered  to  the  shore  station  via  5  to  25 
KHz  bandwidth  FM  transmitters,  utilizing  modified  AN/SSQ- 
57A  sonobuoy  components. 
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The  signals  were  received  at  the  shore  station  by  an  ARR- 
52A  four  channel  sonobuoy  receiver.  Figure  11  schematically 
shows  the  data  acquisition  and  real  time  processing  system 
located  at  the  shore  station.  All  four  hydrophone  signals 
were  fed  to  an  analogue  sample  and  hold  unit,  while  HI  and 
H2  were  also  fed  to  the  dual-channel  spectrum  analyzer. 

Each  channel  of  the  analogue  processor  had  a  quasi¬ 
integrator,  gating  unit,  threshold  detector,  and  sample  and 
hold  system.  Amplitude  and  phase  fluctuations  of  the  20  KHz 
signal  were  simultaneously  measured  by  a  quasi-integration 
of  the  positive  half  of  the  signal.  The  signal  was  passed 
to  the  sample  and  hold  unit  once  a  preset  threshold  was  ex¬ 
ceeded,  during  a  computer  controlled  acquisition  time-window. 
The  sample  window  width  was  set  at  3  ms  in  order  to  minimize 
the  effect  of  ambient  noise  pretriggering  the  sample  and  hold 
unit,  and  thereby  producing  erroneous  data.  Representative 
voltages  were  stored  in  the  sample  and  hold  unit  until  ad¬ 
dressed  by  the  Hewlett-Packard  9845,  where,  one  channel  at 
a  time,  the  signals  were  measured  by  a  scanner  and  five  digit 
DVM,  and  subsequently  logged  into  the  data  acquisition  sys¬ 
tem.  During  the  same  5-second  cycle,  the  atmospheric  data 
(wind  speed  and  direction  and  insolation)  were  also  scanned, 
measured  and  logged  into  the  computer. 

The  5  ms  time  records  for  Hi  and  H2  (the  two  bottom- 
mounted  hydrophones)  were  digitized  to  12  bit  accuracy  by 
the  spectrum  analyzer  after  a  predetermined  digital  delay 
initiated  by  the  pulse  generator  circuitry.  For  each  pulse, 
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DATA  ACQUISITION  AND 
REAL-TIME  PROCESSING  SYSTEM 


acquisition  and  real  time  processing  schematic 


the  two  1024  point,  5  ms  records  were  transferred  to  the 
Kennedy  9832  digital  tape  recorder  via  the  computer.  The 
digital  delay  was  set  at  256.5  ms  (the  travel  time  between 
transmitter  and  receiver  station)  with  an  accuracy  of  100 
Nanoseconds. 

The  schematic  diagram  of  the  real  time  processing  and 
recording  of  oceanic  data  is  also  shown  in  Figure  11.  All 
oceanic  measurements  were  sampled  at  a  rate  of  32  samples 
per  second  by  the  PCM  encoders.  The  remote  stations  (R  &  P) 
were  actuated  every  2  hrs  for  30  min  by  crystal  controlled 
clocks. 

From  each  of  the  offshore  stations,  the  analog  data  was 
sent  to  an  analog  preconditioning  unit  and  then  to  a  PCM 
encoder,  where  a  serial  bi-phase  output,  including  the  data 
and  synchronizing  words,  was  telemetered  to  shore.  The  sig¬ 
nals  were  received  via  a  four  channel  PCM  receiver  system 
(one  channel  per  station)  and  recorded  on  an  eight  channel 
analog  tape  recorder .  At  the  same  time ,  selected  channels 
were  sent  through  a  PCM  decoder  for  real-time  viewing  and 
analysis  on  an  oscilloscope  and  chart  recorder. 

D.  DATA  PROCESSING 

The  four-day  period  from  15  through  18  August  was  chosen 
for  this  preliminary  analysis.  The  acoustic  data  were  digi¬ 
tally  recorded  in  a  format  compatible  with  the  NPS  main¬ 
frame  computer.  However,  it  was  discovered  that  the  Kennedy 
9832  recorder  was  faulty  during  the  experiment,  resulting 


49 


in  some  lost  data  and  other  byte  shifted  data.  The  byte 
shifting  problem,  although  sparse,  required  complicated 
additions  to  the  data  analysis  scheme  to  check  for  byte 
shifts,  and  make  corrections  where  possible.  The  oceanic 
data  was  pre-processed  by  the  HP-9845  computer  to  generate 
IBM  compatible  data  tapes. 

1 .  20  KHz  Acoustics  Data  Processing 

Two  methods  were  used  to  determine  the  amplitude  of 
the  20  KHz  sinewave  pulse  transmitted  at  the  start  of  each 
composite  pulse.  The  first  technique  used  analog  sample  and 
hold  processors  which  measured  in  real-time  the  peak  value 
of  the  positive  half  of  four  cycles  of  the  20  KHz  pulse. 

The  second  estimate  of  the  20  KHz  amplitude  was  obtained  from 
computer  analysis  of  the  digitized  hydrophone  signals. 

During  the  first  stage  of  the  signal  processing,  the 
precise  arrival  time  of  the  pulse  was  determined  from  the 
cross-correlation  between  a  section  of  the  received  pseudo¬ 
random  noise  pulse  and  a  reference  transmitted  signal.  This 
arrival  time  was  used  to  identify  the  beginning  and  end  posi¬ 
tion  of  each  20  KHz  pulse.  The  RMS  energy  of  the  signal  was 
then  calculated  over  8  of  the  20  KHz  cycles  and  then  stored. 
The  time  series  of  acoustic  amplitude  and  travel-time  obtained 
from  these  two  methods  were  plotted  and  stored  on  disk  for 
further  analysis. 

2 .  Oceanic  Data  Processing 

A  data  processing  chain  under  computer  control  was 
developed  to  decode,  deglitch  and  sort  the  PCM  encoded  oceanic 
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data  gathered  at  the  shore,  transmitter  and  receiver  sta¬ 
tions.  Data  from  one  station  at  a  time  was  time  synchronized 
and  the  serial  data  stream  switched  to  the  PCM  decoder.  The 
decoded  data  were  passed  to  a  HP-9845  where  the  channel  iden¬ 
tifiers  were  stripped  and  the  data  deglitched  according  to 
slew-rate  criteria  selected  for  each  channel.  These  data 
were  then  averaged  to  reduce  the  effective  sampling  rate 
from  32  per  second  to  2.7  per  second,  and  the  4096  data  points 
from  each  channel  were  then  temporarily  stored  on  digital 
cassettes.  When  all  three  stations  had  been  preprocessed, 
simple  statistics  were  calculated  and  printed  for  30  selected 
output  channels.  The  statistical  information  and  the  pre- 
processed  data  were  then  transferred  to  IBM-compatible  digital 
tape. 
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IV.  DATA  ANALYSIS 


Four  days,  15-18  August  1979  are  described  and  analyzed 
for  this  preliminary  examination  of  AVEX  results.  The 
weather  during  the  four  days  was  overcast,  with  fog  in  the 
morning  and  only  very  occasional  partial  clearing  (in  the 
afternoon  of  16  August) .  The  winds  were  generally  light, 
with  light  winds  during  the  night  and  morning,  and  in¬ 
creasing  to  10  knots  or  less  in  the  afternoon.  The  baro¬ 
metric  pressure  held  steady  at  30.0  in.  of  mercury. 

The  gross  temperature  profiles  were  measured  with  XBT's 
and  reflected  the  rather  constant  atmospheric  conditions 
during  the  four  days.  The  sea  surface  temperature  varied 
little  from  13 °C.  The  profiles  showed  a  very  gradual  shallow 
thermocline  in  the  upper  10m,  then  a  partially  mired  layer 
down  to  about  50m,  a  thermocline  to  about  100m  and  then 
isothermal  conditions  to  the  bottom  temperature  of  about  10°C. 
The  calculated  Vaisala  frequencies  at  the  30m  depth  of  the 
transmitter  and  receiver,  based  on  the  temperature  profiles, 
were  high,  ranging  between  0.05 — 0.1  hz. 

The  temperatures,  current  speed  and  direction  measured 
by  the  Aanderra  current  meter  at  46m  in  mid-Canyon,  and  tidal 
elevation  for  16  to  18  August  are  shown  in  Figure  12.  The 
tidal  elevations  were  measured  at  Wharf  Number  Two  in  Monterey 
harbor,  10  Km  to  the  north,  and  there  appears  to  be  a  slight 
lag  with  the  current  meter  data.  The  current  speeds  were 
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Figure  12.  Temperature,  current  direction  and  speed,  and  tidal  elevation  measured 
with  the  46m  depth  Aanderra  current  meter  in  Carmel  Canyon. 


generally  weak  and  show  a  definite  tidal  component  with 
considerable  variability  superposed.  The  current  direction 
signal  is  noisy  and  unreliable  because  the  current  speeds 
were  near  the  threshold  for  the  vane.  The  tidal  contribution 
to  the  temperature  signal  is  not  as  obvious,  but  consider¬ 
able  micro-structure  in  the  temperature  signal  can  be 
observed. 

A  log-log,  frequency  spectrum  of  the  temperature  fluctua¬ 
tions  for  the  18-day  period  5-23  August  is  shown  in  Figure  13 
Identifiable  peaks  in  the  spectrum  are  the  inertial  period 
(0.049  cph)  and  the  semidiurnal  tide  (0.081  cph) ,  and  its 
first  and  second  harmonics  (0.161  cph  and  0.243  cph)  At 
higher  frequencies,  the  spectral  level  decreases  in  the 
internal  wave  region.  The  -2  and  -2.5  spectral  slopes  pro¬ 
posed  in  the  Garrett-Munk  (1979)  internal  wave  model  are 
indicated.  At  the  highest  frequencies,  the  spectrum  flattens 
out  into  the  microstructure  region.  The  spectrum  appears 
noisy  because  it  was  calculated  with  only  8  degrees  of  free¬ 
dom  to  resolve  the  inertial  and  tidal  peaks. 

Examples  of  25-min  oceanic  data  starting  on  hours  0100, 
0300,  0700,  and  0900  on  17  August  are  shown  in  Figures  14-17. 
Starting  at  the  top  of  each  figure,  the  measured  variables 
listed  on  the  left  hand  side  of  the  figures  included:  baro¬ 
metric  pressure,  back  radiation,  incident  radiation,  wind 
direction,  wind  speed,  receiver  tower  battery  voltage  monitor 
water  surface  elevation,  thermistors  on  the  receiver  tower 
R12-R1,  ducted  current  meters  on  transmitter  tower  CM2-1, 
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Figure  17.  25-minute  environmental  data  records 
starting  0900  hours,  17  August. 
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transmitter  voltage  monitor,  and  thermistors  mounted  on 
transmitter  tower,  T9-T1.  The  units  and  range  of  plot  values 
are  listed  on  the  righthand  side  of  figures.  Thermistors 
on  the  transmitter  tower,  T9-T1,  have  glitches  due  to  the 
feedback  of  biphase  encoded  signals  into  sensitive  thermistor 
signal  conditioners,  and  were  not  analyzed.  All  the  other 
sensors  shown  in  the  figure  worked  well.  The  temperature 
microstructure  time  series  showed  great  variability.  Focus¬ 
ing  attention  on  the  thermistors  at  the  receiver  tower  in 
Figures  14  to  17,  recall  that  the  vertical  array  consisted 

s. 

of  thermistors  R1  to  5  and  the  horizontal  array  contained 

thermistors  R5  to  12,  with  R5  common  to  both  (see  Figure  7). 

I 

At  0100  hour,  the  mean  temperature  was  relatively  constant 
at  10.5  °C  with  visible  microstructure.  During  0300  hour, 
a  portion  of  an  internal  wave  is  visible  going  from  the 
crest  (colder  water  brought  up)  to  the  trough  (warmer  water) 
with  a  change  of  about  0.2  °C.  The  waves  shapes  are  very 
characteristic  of  nonlinear  internal  waves  propagating  onto 
a  shelf  (see  Phillips,  1977,  p.  199). 

The  time  series  during  0700  hour  appears  to  show  a  short 
period  internal  wave  trough  upon  which  extensive  microstruc¬ 
ture  is  observed.  An  internal  wave  going  from  the  trough 
to  crest  appears  during  the  0900  hour  record  with  very  ex¬ 
tensive  microstructure  and  high  frequency  internal  waves 
occurring.  The  gross  temperature  variations  appear  consistent 
.  with  internal  waves  slowly  propagating  by  the  array,  and 

less  correlated  microstructure,  that  at  times  is  reminiscent 
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of  billow  turbulence  caused  by  wave  breaking.  Propagation 
velocity  can  be  calculated  from  the  thermistor  arrays  by 
observing  the  travel  time  of  the  wavefront  across  the  array. 
Referring  to  Figure  17,  for  example,  a  line  connecting  the 
wavefront  across  the  seven  horizontal  thermistors  indicates 
a  time  of  about  100  seconds  to  travel  the  340cm,  which 
corresponds  to  a  velocity  of  3  to  4  cm/sec;  this  speed  is 
consistent  with  the  velocities  measured  in  the  center  of  the 
Canyon  with  the  Aanderra  current  meter. 

Correlation  lengths  were  measured  from  the  calculated 
autocorrelation  functions.  Examples  of  structure  functions 
and  autocorrelation  functions  calculated  for  the  same  hori¬ 
zontal  and  vertical  thermistor  arrays  at  the  hours  0100,  0300, 
0700,  and  0900  on  17  August  are  shown  in  Figures  18  and  19. 

The  structure  functions  and  autocorrelation  functions  were 
calculated  using  all  combinations  of  thermistors  possible. 

The  heavy  horizontal  bar  indicates  the  length  of  the  arrays. 

In  general,  the  arrays  were  too  short  to  completely  measure 
the  structure  and  autocorrelation  function,  particularly  in 
the  horizontal.  Integral  scale  lengths,  corresesponding  to 
the  average  eddy  size,  were  calculated  from  the  area  under 
the  autocorrelation  functions.  The  integral  scale,  l,  is 
defined 

00 

/  R„(p)dp 

*  =  °“Rx-(6)  (21) 
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Figure  19.  Vertical  array  structure  funct1ons(top  row)  and  autocorrelation 
funct1ons(second  row),  and  horizontal  array  structure  funct1ons( third  row) 
and  autocorrelation  functions (bottom  row)  for  0700( 1  eft  column)  and  0900 
(right  column),  17  August. 
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When  the  array  was  too  short,  for  example  at  0100,  the 
measured  autocorrelation  functions  were  linearly  extrapo¬ 
lated  to  a  zero-crossing.  The  integral  scale  is  a  rela¬ 
tively  rough  measure  and  is  not  very  sensitive  to  the 
extrapolation  procedure. 

The  Taylor  micro-scale  can  also  be  calculated  by 
assuming  the  micro-scale  is  quasi-Gaussian ,  where  the  auto¬ 
correlation  function  is  described  by 

2  2 

Vp)  =  Rx(<l>)e'6  p  (22) 


where  S  is  a  decay  constant. 

The  microscale,  M,  for  a  turbulent  medium  is  defined 
by  fitting  a  parabola  to  the  smallest  scales  and  measuring 
the  zero-cross,  and  is  given  by 


d2Rx(p) 

dx2 


p-*-0 


(23) 


For  a  process  described  by  (23) , 


M  =  1.12  i 


(24) 


Examining  first  the  examples  of  the  vertical  array  autocorre¬ 
lation  functions,  at  0100,  the  integral  scale  length  was 
l  =  570cm;  at  0300,  the  correlation  length  decreased  to 
£  =  152cm,  £  =  230cm  at  0700,  and  £  very  large  at  0900.  The 
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horizontal  array  correlation  lengths  were  always  greater, 
with  i  =  760cm  at  0100,  l  very  large  at  0300,  l  =  1,560cm  at 
0700,  and  i  =  860cm  at  0900.  The  integral  scales  for  the 
vertical  and  horizontal  arrays  for  the  four  days  are  plotted 
later  (Figures  25  and  26) .  The  integral  scales  varied 
considerably  with  time.  The  anisotropy  of  the  vertical  to 
horizontal  scales  is  evident,  but  not  very  strong,  with  the 
average  Z  (horizontal) /£  (vertical)  2.  4. 

Five  examples  of  half-hour  time  series  of  pulse  travel 
time  are  shown  in  Figure  20.  Each  plot  contains  three  hundred 
travel-time  estimates  calculated  from  the  digitized  hydro¬ 
phone  signal,  and  examples  have  been  selected  from  August  16, 
hours  0100,  0300,  0500,  0700  and  1900.  The  travel  time 
plots  are  automatically  scaled  to  show  details  of  the  varia¬ 
tions,  and  have  arbitrary,  but  fixed,  mean  offsets  for  each 
channel.  These  plots  consistently  show  slow  and  relatively 
smooth  changes  in  the  acoustic  travel  time,  with  fluctuation 
periods  of  the  same  order  as  the  half-hour  sample  period. 

In  periods  of  low  activity,  such  as  Figure  20,  the  scatter 
of  the  one  microsecond  measurement  resolution  steps  can  be 
seen,  suggesting  that  the  rms  noise  of  the  travel  time  esti¬ 
mates  is  approximately  one  microsecond.  The  standard  deviation 
of  the  travel  time  fluctuations  were  calculated  for  each  half 
hour  sample  between  August  16  and  18,  and  were  found  to  range 
between  1  and  60us  for  both  hydrophones.  These  standard 
deviation  estimates  are  plotted  in  Figure  21.  It  is  clear 
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Figure  20.  Travel-time  estimates  for  16  August: 

a )  Hour  01 00 

b)  Hour  0300 

c)  Hour  0500 

d)  Hour  0700 

e)  Hour  1 900 


from  the  scatter  of  the  points  that  there  are  large  changes 
in  the  variability  of  sound-speed  propagation  across  the 
Canyon  between  each  half-hour  sampling  period.  Some  of  this 
scatter  in  the  standard  deviation  estimates  is  due  to  under¬ 
sampling  of  partially  periodic  travel  time  fluctuations,  but 
inspection  of  other  time-series  in  the  four-day  sample  show 
long  quiescent  periods  followed  by  bursts  of  high  sound- 
speed  variations. 

The  corresponding  20  KHz  pulse  amplitude  time-series  are 
shown  in  Figure  22.  The  most  significant  differences  between 
the  travel-time  and  amplitude  data  are  more  rapid  fluctuations 
in  the  amplitude  signal.  Additionally,  there  were  signifi¬ 
cant  differences  in  the  amplitude  received  by  the  two  hydro¬ 
phones,  which  were  separated  by  only  2m;  this  suggests  that 
the  intensity  fluctuations  are  caused  by  refractive  index 
fluctuations  with  scales  much  shorter  than  those  associated 
with  integral  sound-speed  variations  recorded  in  the  travel 
time  measurement.  Flatte  (1979)  has  suggested  that  small 
changes  of  the  geometry  in  the  sound  speed  fluctuations  along 
the  acoustic  ray  path  cause  focusing  and  defocusing  of  the 
rays  which  are  seen  as  large  intensity  variations  at  the  re¬ 
ceiver.  The  measurements  of  rapidly  changing  intensity  made 
in  this  experiment  support  this  hypothesis. 

To  examine  the  relationship  between  travel  time  and 
amplitude  variations,  a  correlation  plot  was  constructed 
(Figure  23) .  Despite  the  large  scatter,  the  correlation 


63 


r:ie  \minj 


;*e  mm 


Figure  22.  20KHz  pulse  amplitude  for  16  August: 

a)  Hour  0100 

b)  Hour  0300 

c)  Hour  0500 

d)  Hour  0700 

e)  Hour  1900 


plot  does  suggest  that,  while  large  intensity  fluctuations 
occurred  without  large  integral  sound-speed  variations,  the 
converse  is  not  true,  i.e.,  large  intensity  fluctuations 
occurred  in  periods  with  minimal  travel  time  changes. 

To  place  the  AVEX  results  in  Flatte's  strength/ 
diffraction  space  (see  equations  18  and  20)  the  values  of 
$  and  A  were  calculated  from  the  correlation  lengths  measured 
at  the  receiver  array  for  an  acoustic  frequency  of  20  KHz. 
These  estimates  are  plotted  in  Figures  24  and  25,  using 
correlation  lengths  measured  at  the  vertical  and  horizontal 
thermistor  arrays,  respectively.  Although  equations  18  and 
20  do  not  strictly  apply  in  the  anisotropic  conditions  found 
at  the  AVEX  site,  the  very  low  propagation  angle  (due  to 
source  and  receiver  being  at  the  same  depth  ±2m)  allows  us 
to  assume  isotropic  conditions  (see  e.g.,  Desaubie,  1979). 

Even  though  the  horizontal  and  vertical  scales  differ  signi¬ 
ficantly,  the  points  plotted  in  Figures  24  and  25  show  the 
approximate  position  in  the  <t>-A  space  of  each  half  hour  acous¬ 
tic  measurement.  With  two  exceptions,  all  the  points  lie 
well  within  the  geometrical  optics  regime. 

The  relationship  between  the  temperature  lengths  scales, 
temperature  variance,  and  the  acoustic  pulse  amplitude  may 
be  seen  in  the  time-series  plotted  in  Figure  26.  These  plots 
fail  to  suggest  any  obvious  relationship  between  the  varia¬ 
bles,  but  does  show  the  large  aperiodic  fluctuations  over 
the  sixty-hour  sample  interval.  A  similar  set  of  time-series 
including  travel-time  standard  deviations  is  shown  in  Figure 
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Figure  24.  Strength,  *,  versus  diffraction.  A,  for  16-18  August, 
using  vertical  length  scales  measured  at  the  receiver 
tri pod . 
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27.  Again  the  large  variations  of  the  travel-time  can  be 
seen. 

To  test  whether  a  causal  relationship  exists  between  the 

measured  temperature  scales  and  variance,  correlation  plots 

using  Chernov's  criteria  and  the  Debye  Approximation  were 

made.  In  Figure  28,  pulse  amplitude  is  plotted  against 
-3/2 

where  aT  is  the  temperature  standard  deviation  and 
is  the  horizontal  integral  scale  length,  to  test  Chernov's 
formula  for  amplitude  variance  (equation  12) .  The  large 
scatter  in  this  plot  indicates  that  the  relationship  is  not 
confirmed  for  these  measurements. 

A  similarly  poor  correlation  was  obtained  when  amplitude 
variation  was  plotted  against  a  /XT  (see  Figure  29)  to  test 
the  Debye  approximation  (equation  16) . 

The  reason  for  the  very  large  scatter  in  these  correla¬ 
tion  plots  is  probably  that  both  the  length-scales  and 
temperature  variances  measured  at  the  receiver  tower  were 
not  always  typical  of  the  temperature  field  along  the  entire 
transmission  path.  (The  failure  of  the  yo-yo  profiler,  which 
was  to  have  provided  the  principle  temperature  structure 
data  at  about  mid-range,  prevents  better  estimates  of  the 
refractive  index  field.) 
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a)  Horizontal  Integral  scale(m)  vs.  time(hours), 

b)  Temperature  standard  deviat1on(°C )  vs.  tlme(hours), 

c)  Vertical  length-scale(m)  vs.  time  (all  measured  at 
the  receiver  array) , 

d)  Travel-time  standard  deviation^ s) . 
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Figure  28.  Correlation  between  amplitude  variations,  length-scales, 
and  temperature  variance  using  Chernov's  model  (arbitrary 
units) . 
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Figure  29.  Correlation  between  amplitude  variations,  length-scales, 
and  temperature  variance  using  the  Debye  approximation 
(arbitrary  units). 
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V.  SUMMARY  AND  CONCLUSIONS 


The  Acoustic  Variability  Experiment  (AVEX)  described 
here  measured  acoustic  pulse  propagation  along  a  wholly  re¬ 
fracted  400m  path  between  a  fixed  source  and  receiver  in  the 
upper  ocean.  The  acoustic  pulse  was  composed  of  a  0.5ms 
burst  of  20  KHz  CW  followed  by  a  4.5ms  pseudo  random  noise 
pulse  containing  energy  in  the  frequency  range  between  5 
and  25  KHz.  The  source  and  receivers  were  stably  mounted 
on  towers  located  on  each  side  of  Carmel  Canyon  at  a  depth 
of  30m. 

The  temporal  variability  of  acoustic  amplitude  and  phase 
was  measured  and  a  preliminary  analysis  of  the  20  KHz  pulse 
data  has  been  described.  The  travel  time  across  the  Canyon 
had  rms  variations  between  1  and  60ys,  calculated  over  each 
half-hour  sample  period,  and  corresponding  20  KHz  pulse  ampli¬ 
tude  fluctuations  between  3%  and  14%.  The  level  of  these 
fluctuations  changed  dramatically  in  the  1.5  hours  between 
successive  measurements. 

Both  the  temporal  and  spatial  structure  of  the  tempera¬ 
ture  field  were  measured  simultaneously  with  the  acoustic 
transmission  at  the  source,  receiver  and  at  a  mid-point. 

A  spectrum  of  a  15-day  temperature  record  showed  the  large 
scale  variability  due  to  inertial  and  tidal  motions  and  to 
internal  waves  at  higher  frequencies.  The  microstructure 
was  intermittent  and  appeared  to  be  associated  with  the 
internal  waves.  The  temperature  integral  scales  calculated 


from  the  spatial  autocorrelation  functions  ranged  from  1.2 
to  14m  (est)  in  the  vertical  and  from  3.4m  to  greater  than 
50m  (est)  in  the  horizontal.  The  temperature  structure  was 
anisotropic  with  the  average  horizontal  to  vertical  scale 
ratio  equal  to  about  four. 

An  attempt  to  compare  the  measured  acoustic  amplitude 
variance  with  temperature  scale  lengths  and  variance  using 
the  theoretical  predictions  of  Chernov  and  Debye  failed  to 
show  any  consistent  relationship.  These  tests  showed  that 
poor  correlations  existed  between  acoustic  intensity  and  the 
temperature  structure  measured  at  the  side  of  the  Canyon  by 
the  receiver  thermistor  arrays.  Due  to  the  failure  of  the 
yo-yo  profiler,  more  representative  measurements  of  the 
temperature  structure  along  the  propagation  path  were  not 
available.  The  experiment  did  demonstrate  the  highly  epi¬ 
sodic  character  of  the  f inestructure  in  the  Canyon  tempera¬ 
ture  field,  which  appeared  to  be  dominated  by  internal 
wave  activity.  The  intermittent  bursts  of  activity  were 
reflected  in  both  the  thermistor  array  temperature  measure¬ 
ments,  and  indirectly  by  the  pulse  intensity  and  travel-time 
estimates. 
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